Scattering theory of plasmon-assisted entanglement transfer and distillation 
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We analyse the quantum mechanical limits to the plasmon-assisted entanglement transfer observed 
by E. Altewischer, M.P. van Exter, and J.P. Woerdman [Nature, 418, 304 (2002)]. The maximal 
violation S of Bell's inequality at the photodetectors behind two linear media (such as the perforated 
metal films in the experiment) can be described by two ratio's n, T2 of polarization-dependent 
transmission probabilities. A fully entangled incident state is transferred without degradation for 
T\ = T2, but a relatively large mismatch of n and T2 can be tolerated with a small reduction of S. 
We predict that fully entangled Bell pairs can be distilled out of partially entangled radiation if n 
and T2 satisfy a pair of inequalities. 
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The motivation for this work came from the recent re- 
markable demonstration by Altewischer, Van Exter, and 
Woerdman of the transfer of quantum mechanical entan- 
glement from photons to surface plasmons and back to 
photons . Since entanglement is a highly fragile prop- 
erty of a two-photon state, it came as a surprise that this 
property could survive with little degradation the con- 
version to and from the macroscopic degrees of freedom 
in a metal Q. 

We present a quantitative description of the finding 
of Ref. [| that the entanglement is lost if it is mea- 
sured during transfer, that is to say, if the medium 
through which the pair of polarization-entangled photons 
is passed acts as a "which-way" detector for polarization. 
Our analysis explains why a few percent degradation of 
entanglement could be realized without requiring a highly 
symmetric medium. We predict that the experimental 
setup of Ref. Q could be used to "distill" % || fully 
entangled Bell pairs out of partially entangled incident 
radiation, and we identify the region in parameter space 
where this distillation is possible. 

We assume that the medium is linear, so that its ef- 
fect on the radiation can be described by a scattering 
matrix. The assumption of linearity of the interaction of 
radiation with surfaceplasmons is central to the litera- 
ture on this tòpic 0, la, 0, H m ■ We will not make any 
specihc assumptions on the mode and frequeney depen- 
dence of the scattering matrix, but extract the smallest 
number of independcntly mcasurable parameters needed 
to describe the experiment. By concentrating on model- 
independent results we can isolatc the fundamental quan- 
tum mechanical limitations on the entanglement transfer, 
from the limitations specific for any particular transfer 
mechanism. 

The system considered is shown schematically in Fig. 
^ Polarization-entangled radiation is scattered by two 
objeets and detected by a pair of detectors behind the 
objeets in the far-ficld. The objeets used in Ref. |l| 
are metal films perforated by a square array of subwave- 
length holes. The transmission amplitude t aa i,i of object 
i = 1,2 relates the transmitted radiation (with polar- 



ization a = H, V) to the incident radiation (polarization 
a' = H, V). We assume a single-mode incident beam and 
a single-mode detector (smaller than the coherence area), 
so that we require a set of eight transmission amplitudes 
t Ga i i out of the entire scattering matrix (which also con- 
tains reflection amplitudes and transmission amplitudes 
to other modes). The extension to a multi- mode theory 
(needed to describe some aspeets of the experiment 0) 
is left for a future investigation. 
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FIG. 1: Main plot: Efficiency of the entanglement transfer 
for a fully entangled incident state, as given by Eq. 114^ . 
The maximal violation Smax of Bell's inequality at the pho- 
todetectors is plotted as a function of the ratio ri/ra = 
T1+T2- /T1-T2+ of the polarization-dependent transmission 
probabilities. The inset shows schematically the geometry 
of the experiment 0. A pair of polarization-entangled pho- 
tons is incident from the left on two perforated metal films. 
The photodetectors at the right, connected by a coincidence 
counter, measure the degree of entanglement of the transmit- 
ted radiation. 



2 



We do not require that the scattering matrix is uni- 
tary, so our results remain vàlid if the objects absorb 
part of the incident radiation. What is neglected is the 
thermal radiation, cither from the two objects or from 
the electromagnètic environment of the detectors. This 
thermal noise is insignificant at room temperature and 
optical freqüències. 

The radiation incident on the two objects is in a known, 
partially entangled state and we wish to determine the 
degree of entanglement of the detected radiation. It is 
convenient to use a matrix notation. The incident two- 
photon state has the general form 

|* in > = a S H |HH)+al 1 n v |HV)+ a V n H|VH)+ a ^ v |VV). (1) 

The four complex numbers , form a matrix 



-4; 



"HH a HV 

in in 
"VH "VV 



(2) 



Normalization of l^in) requires 
the trace of a matrix. 

The four transmission amplitudes í CT0 -',i of object i 
1 2 form the matrix 



Ti = 



(3) 



The transmitted two-photon state |^ ut) has matrix of 
cocfficients 



A out = Z-^T^Tl 
with normalization factor 

Z = Tr(T 1 A in T^)(T 1 A in T^. 



(4) 



(5) 



(The superseript "t" denotes the transpose of a matrix.) 

We quantify the degree of entanglement in terms of the 
Clauser-Horne-Shimony-Holt parameter S |10|. which 
measures the maximum violation of Bell's inequality and 
was used in the experiment pj. This parameter can be 
obtained from a decomposition of into a superposi- 
tion of a fully entangled state (with weight y/P) and a 
factorized state orthogonal to it 0, 0] . The relation is 



S = 2\Jl + P 2 , P 2 =4DetAA t , 



(6) 



with "Det" the determinant and < P < 1. (The con- 
currence 01 * s identical to P.) A fully entangled state 
has P = 1, S = 2v2, while a factorized state has P = 0, 
S = 2. The fully entangled state could be the Bell pair 
(|HV) — |VH))/\/2~, or any state derived from it by a 
local unitary transformation (A — » U AV with £/, V ar- 
bitrary unitary matrices). The degree of entanglement 
p n = 2|DetAi n | of the incident state is given and we 
seek the degree of entanglement P ou t = 2|Det A out \ of the 
transmitted state. We are particularly interested in the 
largest P ou t that can be reached by applying local uni- 
tary transformations to the incident state. This would 



correspond to the experimental situation that the polar- 
izations of the two incoming photons are rotated inde- 
pendently, in order to maximize the violation of BclPs 
inequality of the detected photon pair. 

Before proceeding with the calculation we introduce 
some parametrizations. The Hermitian matrix product 
T^P/ has the eigenvalue-eigenvector decomposition 



T x Tl 



C/T 



Ti 



i+ 







Ti_ 



u, r a rj 



2+ 







To 



V. 

The matrices of eigenvectors U, V are unitary and the 
transmission eigenvalues Tj± are real numbers between 
and 1. We order them such that < X! t _ < T i+ < 1 
for each i = 1, 2. We will see that the maximal entangle- 
ment transfer depends only on the ratios Tj = T; + /Ti_. 
This parametrization therefore extracts the two signifi- 
cant real numbers t\ , r 2 out of cight complex transmis- 
sion amplitudes. The Hermitian matrix product ^4; n ^4; n 
has eigenvalues A± = i± i(l — P^) 1 ^ 2 - These appear in 
the polar decomposition 



UA in V 



U+ U- 

—u*_ u*. 



X A 4 








V-\- V- 



(8) 

The phase (f> is real and u±,v± are complex numbers con- 
strained by |w±| = (\±u) 1 / 2 , \v±\ = (| ±v) 1 / 2 , with real 
u, v e (— i, i). These numbers can be varied by local uni- 
tary transformations, so later on we will want to choosc 
vàlues which maximize the detected entanglement. 

With these parametrizations a calculation of the de- 
terminant of A out leads to the following relation between 
P in and Pout: 



Pont — 



PnV T l T 2 



(ri-l)(r 2 -l) 



-[A+Q+ + A_Q_ 



2v /A^(± - U Y /2 (| " « 2 ) 1/2 cos$r\ (9) 



Q± 



u± 



ri + 1 
] ti-1 



v± 



1 T 2 



T2 



1 



(10) 



The phase $ equals the argument of u + u*_v + v- . To max- 
imize P out we should choose $ = 0. 

We first analyze this expression for the case of a fully 
entangled incident state, as in the experiment Q. For 
p n = 1 one has A + = A_ = 1/2, and Eq. (JÜJl simplifies 
to 

P 

(T 1 + l)(T 2 + l) + 4a(r 1 -l)(r 2 -l)' 1 ' 
a = uv-{\- u 2 ) x / 2 {\ - v 2 ) 1 ' 2 cos $. (12) 

Since Ti > 1 and \a\ < | we conclude that the degree of 
entanglement is bounded by P m i n < Pout < Pnax, with 



Pn 



1 + TlT 2 ' 



Pn 



2^7 



T-2 



i + n/T 2 ' 



(13) 
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The maximum P max can always be reached by a proper 
choice of the (fully entangled) incident state, so the max- 
imal violation of Bell's inequality is given by 



<? — 9 




(14) 



The dependence of S mSíX on T\jvx is plotted in Fig. ^ 
Full cntanglement is obtained for n = ts, hence for 
Ti + T 2 _ = Ti_T 2+ . Generically, this requires either iden- 
tical objects (Ti± = Ï2±) or non-identical objects with 
Xí + = Tj_. If Ti = T2 there are no "which-way" labels 
and entanglement fully survives with no degradation. 

Small deviations of T1/T2 from unity only reduce the 
cntanglement to second order, 



5 max = 2V2[1 - 4( T i/ r 2 - l ? 
lo 



+ 0(r 1 /r 2 -l) 3 l. (15) 



So for a small reduction of the entanglement one can 
tolerate a large mismatch of the transmission probabil- 
ities. In particular, the experimental result S = 2.71 
for plasmon-assisted entanglement transfer .1] can be 
reached with more than a factor two of mismatch (S = 
2.71 for ti/t 2 = 2.4). 

As a simple example we calculate the symmetry pa- 
rameter ti/t 2 for a Lorentzian transmission probability, 
approjjriate_ for plasmon-assisted entanglement transfer 
We take 



T, 



TV 2 



(üj - lü í± ) 2 + r 2 ' 



(16) 



where lüq is the frequeney of the incident radiation, T 
is the linewidth, and T is the transmission probabil- 
ity at the resonance frequeney u>i±. (For simplicity we 
take polarization-independent T and T.) The transmis- 
sion is through an optically thick metal film with a rect- 
angular array of subwavelength holcs (lattice constants 
Li±). The dispersion relation of the surface plasmons is 
lüí± = (1 + l/e) 1 / 2 2-!rnc/Li± where e is the real part 
of the dielectric constant and n is the order of the reso- 
nance, equal to the number of plasmon-field oscillations 
in a lattice constant. We break the symmetry by taking 
one square array of holes and one rectangular array (lat- 
tice constants Lq = L\+ = L 2 + = £2- and L\ = L%-). 
The lattice constant Lç, is chosen such that the incident 
radiation is at resonance. The symmetry parameter be- 
comes 



T2 



= 1 + (27T) 5 



nl 

To 



nl 
LÍ 



l = 



1 



(17) 



The length l is the propagation length of the surface plas- 
mon. (We have taken c(l + l/e) 1 / 2 for the plasmon group 
velocity, vàlid if lüq is not close to the plasma frequeney 
9].) Combining Eqs. (O and O we see that the 
deviation of S'max from 2y/2 (the degradation of the en- 
tanglement) is proportional to the fourth power of the 




FIG. 2: The shaded strips indicate the vàlues of ln ri and ln T2 
for which P out = 1 can be reached with P[ n = 0.5 (horizontally 
shaded) and Pi n = 0.9 (vertically shaded), in accordance with 
Eq. fljp . 



difference between the number of oscillations of the plas- 
mon field along the two lattice vectors. 

Turning now to the more general case of a partially 
entangled incident state, we ask the following question: 
Is it possible to achieve P ou t = 1 even if P\ n < 1? In 
other words, can one detect a 2^2 violation of Bell's 
inequality after transmission even if the original state 
was only partially entangled? Examination of Eq. (JÜJ) 
shows that the answer to this question is: Yes, provided 
T\ and t 2 satisfy 



ln — I < 2 arcosh(P ir 
t 2 



L ) and ln nr 2 > 2 arcosh(P in 1 ). 

(18) 

The allowed vàlues of t\ and r 2 lie in a strip that is open 
at one end, see Fig. [3 The boundaries are reached at 
\u\ = \v\ = h. The region inside the strip is reached by 
choosing both |it| and |u| < 1/2. For P in = 1 the strip 
collapses to the single line t\ = r 2 , in agreement with Eq. 

03- 

The possibility to achieve P out = 1 for Pn < 1 is 
an example of distillation of entanglement Hl|. (See 
Ref. [lA ha lï^ . IÏ7l for other schemes proposed re- 
cently, and Ref. [Ï8J for an experimental realization.) 
As it should, no entanglement is created in this oper- 
ation. Out of iV incoming photon-pairs with entangle- 
ment p n one deteets NZ pairs with entanglement P ou t = 
flnZ-Vïi+Ti-Ta+Ta-, so that NZP out < NP in . 

In conclusion, we have shown that optical entangle- 
ment transfer and distillation through a pair of lin- 
ear media can be described by two ratios T\ and t 2 of 
polarization-dependent transmission probabilities. For 
fully entangled incident radiation, the maximal violation 
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of Bell's inequality at the detectors is given by a func- 
tion l|14|) of ti/t 2 which decays only slowly around the 
optimal value t\/tz — 1. Distillation of a fully entangled 
Bell pair out of partially entangled incident radiation is 
possible no matter how low the initial entanglement, pro- 
vided that t± and t 2 satisfy the two inequalities Ijl8|l . 

Our results provide a simple way to describe the ex- 
periment Q on plasmon-assisted entanglement transfer, 
in terms of two separately measurable parameters. By 
changing the square array of holes used in Ref. into a 
rectangular array (or, equivalently, by tilting the square 
array relative to the incident beam), one can move away 



from the point n = Ti = 1 and search for the en- 
tanglement distillation predicted here. The possibility 
to extract Bell pairs by manipulating surface plasmons 
may have interesting applications in quantum informa- 
tion processing. 

This work was supported by the Dutch Science Foun- 
dation NWO/FOM and by the US Army Research Of- 
fice (grant DAAD 19-02-0086). J.T. acknowledges the fi- 
nancial support provided through the European Commu- 
nity's Human Potential Program under contract HPRN- 
CT-2000-00144, Nanoscale Dynamics. We have benefit- 
tcd from discussions with J. Preskill and J.P. Woerdman. 



[1] E. Altewischer, M.P. van Exter, and J.P. Woerdman, Na- 

ture 418, 304 (2002). 
[2] W. Barnes, Nature 418, 281 (2002). 
[3] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schu- 

macher, Phys. Rev. A 53, 2046 (1996). 
[4] M.A. Nielsen and I.L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

2000). 

[5] U. Schroter and D. Heitmann, Phys. Rev. B 58, 15419 
(1998). 

[6] M.M.J. Treacy, Appl. Phys. Lett. 75, 606 (1999). 

[7] J.A. Porto, F.J. García-Vidal, and J.B. Pendry, Phys. 

Rev. Lett. 83, 2845 (1999). 
[8] L. Martín-Moreno, F.J. García-Vidal, H.J. Lezec, K.M. 

Pellerin, T. Thio, J.B. Pendry, and T.W. Ebbesen, Phys. 

Rev. Lett. 86, 1114 (2001). 
[9] H. Raether, Surface Plasmons (Springer, Berlin, 1988). 



[10] J.F. Clauser, M.A. Horne, A. Shimony, and R.A. Holt, 

Phys. Rev. Lett. 23, 880 (1969). 
[11] A.F. Abouraddy, B.E.A. Saleh, A.V. Sergienko, and M.C. 

Teich, Phys. Rev. A 64, 050101 (2001). 
[12] N. Gisin, Phys. Lett. A 154, 201 (1991). 
[13] S. Hill and W.K. Wootters, Phys. Rev. Lett. 78, 5022 

(1997); W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[14] H.-W. Lee and J. Kim, Phys. Rev. A 63, 012305 (2000). 
[15] D. Gottesman and J. Preskill, Phys. Rev. A 63, 022309 

(2001). 

[16] T. Yamamoto, M. Koashi, and N. Imoto, Phys. Rev. A 
64, 012304 (2001). 

[17] J.-W. Pan, C. Simón, C. Brukner, and A. Zeilinger, Na- 
ture 410, 1067 (2001). 

[18] P.G. Kwiat, S. Barraza-Lopez, A. Stefanov, and N. Gisin, 
Nature 409, 1014 (2001). 



